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Abstract 

The computation of direct CP asymmetries in charmless B decays at next-to-next- 
to-leading order (NNLO) in QCD is of interest to ascertain the short-distance con¬ 
tribution. Here we compute the two-loop penguin contractions of the current- 
current operators Qi ^2 and provide a first estimate of NNLO CP asymmetries in 
penguin-dominated h ^ s transitions. 


1 Introduction 


Non-leptonic exclusive decays of B mesons play a crucial role in studying the CKM mech¬ 
anism of quark flavour mixing and in quantifying the phenomenon of CP violation. Direct 
CP violation is related to the rate difference of 5 —?■ / decay and its CP-conjugate and 
arises if the decay amplitude is composed of at least two partial amplitudes with different 
re-scattering (“strong”) phases, which are multiplied by different CKM matrix elements. 
Very often useful information on the CKM parameters including the CP-violating phase 
can be obtained from combining different decay modes, whose partial amplitudes are 
related by the approximate flavour symmetries of the strong interaction [1], which are 
then determined from data. 

The direct computation of the partial amplitudes is a complicated strong interaction 
problem, which can, however, be addressed in the heavy-quark limit. The QCD fac¬ 
torization approach [H] employs soft-collinear factorization in this limit to express the 
hadronic matrix elements in terms of form factors and convolutions of perturbative ob¬ 
jects (hard-scattering kernels) with non-perturbative light-cone distribution amplitudes 
(LCDAs). At leading order in A/rrib, 

{M^M2mB) = tml Tl{u) + {M, ^ M,) 


+ 


du / dudv Ti (cu, u, u) fMifpMi(v) (u) 


( 1 ) 


where Qt is a generic operator from the effective weak Hamiltonian. At this order the 
re-scattering phases are generated at the scale nib only, and reside in the loop corrections 
to the hard-scattering kernels. Beyond the leading order factorization does not hold, 
and re-scattering occurs at all scales. The leading contributions to the strong phases are 
therefore of order as{nib) or/and A/nib. It is of paramount importance for the predictivity 
of the approach for the direct CP asymmetries to know whether the short-distance or 
long-distance contribution dominates in practice, since apart from being parametrically 
small, both could be numerically of similar size. 

The short-distance contribution to the direct CP asymmetries is fully known only 
to the first non-vanishing order (that is, 0{as)) through the one-loop computations of 
the vertex kernels T/ performed long ago [21I11IS]. A reliable result presumably requires 
the next-to-next-to-leading order 0(a'^) hard-scattering kernels, at least their imaginary 
parts. For the spectator-scattering kernels T/'^ this task is already completed, both for 
the tree [6H8] and penguin ISlllo] amplitudes, but for the so-called form factor term(s) in 
the first line of ([1]) an important piece is still missing, which is the focus of this Letter. 

We recall that due to CKM unitarity, the amplitude for a B decay governed by the 
b ^ D {D = d,s) transition can always be written in the form 


A(B ^ /) = Af > [T + ...] + Af> [P, + ...], (2) 

where Ap^ — It is generic that the first CKM structure is dominated by the 

colour-allowed or colour-suppressed topological tree amplitude, both denoted by T here. 
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corresponding to the flavonr qnantnm nnmbers of a 6 —)■ uuD transition, while the 
second is dominated by the topological QCD pengnin amplitnde of the b —)■ J2q=u d s 
transition. The hrst is typically larger than the second for D = d and vice-versa for the 
D = s case, which therefore refers to the pengnin-dominated decays snch as B ^ irK 
and related. In the notation of mu, Pc corresponds to the qnantity a 4 (MiM 2 )EI 

The vertex kernels T/ have been compnted at the two-loop 0{a‘^) order only for the 
topological tree amplitndes T [IMS]. However, direct CP asymmetries can only be 
non-zero dne to the interference of the two terms in Q, hence the pengnin amplitnde Pc 
is also needed. Only the one-loop 0{ag) contribntion from the chromomagnetic dipole 
operator Qsg to Pc has been considered in the past na, while the dominant, gennine two- 
loop contribntions remain to be compnted. This calcnlation is technically very involved 
since it reqnires the compntation of massive two-loop pengnin diagrams - a gennine two- 
loop, two-scale problem. One step towards this goal was recently achieved in [15], where 
analytic resnlts of all occnrring master integrals were derived. 

At this point it is important to note that the topological tree and pengnin amplitndes 
are not in one-to-one correspondence with the tree (or cnrrent-cnrrent) operators Q 12 
and QCD pengnin operators Qs-q of the weak effective Hamiltonian. By contracting 
the pp helds of the operators Q ^2 (0 below), they contribnte to the QCD pengnin 

amplitnde starting from the one-loop order. Since these “pengnin contractions” of the 
cnrrent-cnrrent operators come with the large short-distance coefficients Ci^ 2 , we may 
argne that they constitnte the largest contribntion to the pengnin amplitnde at any given 
loop order@ At next-to-leading order we hnd for the pengnin contractions (inclnding the 
chromomagnetic dipole operator Q^g) 

= (-0.0087 - 0.0172^)1^^ ^ + (0.0042 + 0.004H)|^^_^ + 0.0083,^^^, 

a^(7rJ?)|,,o = (-0.0131 - 0.0102^),^^ ^ + (0.0042 + 0.004H)|^^_^ + 0.0083,^^^, (3) 

where we separated the contribntions from the cnrrent-cnrrent and the other operators. 
While there is a cancellation for the real part, the imaginary part from Q 42 is clearly 
dominant. If we add the vertex contractions at leading (LO) and next-to-leading order 
(NLO) and consider the entire form factor contribntion to a\{MiM 2 ) at NLO, the second 
term changes to (—0.0266 -|- 0.0032 Q|q^ ^ in both expressions, and the imaginary part 
from Q ^2 is still by far dominant. Thns, at NLO, the short-distance direct CP asymme¬ 
tries are mainly determined by the one-loop pengnin contractions of the cnrrent-cnrrent 
operators. It is reasonable to assnme that the insertion of QW at two loops also cap- 
tnres the bnlk of the yet nnknown NNLO form factor contribntion T/ to the pengnin 
amplitndes In this Letter we report the resnlt of this compntation together with 

^a\[MiM2) refers to a generically sub-leading penguin contribution to the term multiplied by the 
CKM factor We also note that 04(MiM2) consists of a leading-power term 04 and a power- 

suppressed term Og [Sj. The calculation reported here concerns the leading-power contribution 04. 

^Since the contribution from ^ alone is not renormalization-group invariant, this statement cannot 
be true in arbitrary schemes nor at arbitrary renormalization scales. What we mean is that the statement 
holds in the conventional MS scheme and with a reasonable choice 0 {mi,) of scale. 
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Figure 1: Two-loop penguin diagrams that contribute to the insertion of the operators 
Qi ’2 (black square). The grey hlled circle denotes the one-loop gluon self-energy. 


a few phenomenological implications. We shall provide more technical details together 
with the remaining contributions from the QCD penguin operators Qs-e, which require 
additional work on infrared subtractions not present for Q 12 , Si future publication. 

2 Outline of the calculation 

The effective weak Hamiltonian for b ^ D transitions {D = d, s) is given by 

^ E KoVvoiciQI + + ...) + h.c.. (4) 

^ p=u,c ^ ^ 

Here and in the following we give explicitly only the definitions pertinent to the current- 
current operators relevant to our calculation. We adopt the Chetyrkin-Misiak-Miinz 
(CMM) operator basis [IB], where the current-current operators are dehned as 

Qi = {PlI^THl) Ql = {piYbL) (5l7/4Pl), (5) 

in terms of left-chiral quark fields = |(1 — 75 ) 0 ’- In dimensional regularization the 
operator basis has to be supplemented by evanescent (vanishing m D = A dimensions) 
operators, for which we adopt the convention of mi- 

At the two-loop level about 70 diagrams contribute to the QCD penguin amplitude, 
but only a subset of two dozens (shown in Fig. [T]) are non-vanishing for the insertion of 
the current-current operators The quark in the fermion loop can either be massless 
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(for p = u) or massive (for p = c). In the massless case the problem involves one non¬ 
trivial scale, the momentum fraction u = 1 — u oi the anti-quark in meson M 2 , and 
the structure is similar to the NNLO calculation of the tree amplitudes [nHI3]. In the 
massive case, however, we are dealing with a genuine two-loop, two-scale problem since 
the hard-scattering kernels depend in addition on the mass ratio Sc = As we 

have already elaborated extensively on the kinematics in [15], we shall not repeat those 
formulae here. 

The calculation is performed in dimensional regularization with D = 4 — 2e, where 
ultraviolet (UV) and infrared (soft and collinear) divergences manifest themselves as poles 
in e. The CMM basis ensures that the NDR scheme with a fully anti-commuting 75 can 
be adopted. The amplitude of the diagrams is then computed using standard multi-loop 
techniques. After a Passarino-Veltman [T2] decomposition of the tensor integrals, the 
scalar integrals are reduced to a small set of master integrals by means of integration- 
by-parts techniques naEU] and the Laporta algorithm [2I1122]. To this end, we use the 
program AIR [23] as well as an in-house routine. 

For the massless up-type operator insertions, the diagrams can be expressed in terms 
of the master integrals that appeared in our former calculations [nHI3|. For the massive 
charm-type insertions, on the other hand, we hnd 29 new master integrals. The compu¬ 
tation of the master integrals constitutes the main technical challenge of the calculation. 
Analytic results for all master integrals have recently been derived in [15], based on a dif¬ 
ferential equation approach in a canonical basis [23]. The canonical basis, together with 
suitably chosen kinematic variables, also catalyses the convolution of the hard-scattering 
kernels with the LCDA. 

After the computation of the bare QCD two-loop amplitude, the hard-scattering ker¬ 
nels are extracted from a matching calculation onto soft-collinear effective theory (SCET). 
The main conceptual challenge in this context is the consistent treatment of evanescent 
and Fierz-equivalent operators in SCET, for which we follow the method employed in |13j . 
The SCET operators have the flavour structure J2qiXDXq){^qhv) where x and ^ denote 
collinear light-quark helds moving in opposite directions. The two-loop diagrams relevant 
to the penguin amplitude are all of the “wrong-insertion” type (see [I3]) and hence 
lead to operators where the fermion indices are contracted in the form J2qi^qXq){XDhv)- 
In Fig. [T] the {^qXq) fermion lines correspond to the solid line on the right side of the 
diagram. In the following we omit the sum over q and the flavour labels on the helds. 
In the CMM basis the fermion line that corresponds to (.^y) carries no 75 matrix. This 
suggests that we use the following basis for the SCET operators: 

tL — — 

Oi = X— (1 -Xh)x^i+0--lb)hv, 

On = • • • 7f ""'x x(l + 75)7±«7±/.2n-27±/.2n-3 • • • (6) 

where we need n up to 2 (strings with three 7 matrices in each bilinear). The operator 
Oi is the only physical SCET operator. It is the same as in [TH], whereas the differ 
by the absence of the 1 — 75 factor to the left of x- The operators On are evanescent 
for n > 1. Oi is Fierz-equivalent to Oi/2 in four dimensions, so we add Oi — Oij2 as 
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another evanescent operator. We also recall that the SCET operators are non-local on 
the light-cone [T3] . 

After operator matching the hard-scattering kernels follow from the bare QCD ampli¬ 
tudes plus subtraction terms from UV counterterms of the operators Qi and the SCET 
operators. The master formulae at LO, NLO, and NNLO read, respectively, 


= A} 


(0) 




7 ^( 1 ) _ 4(1)11^ 7(1) 4 (*^) 

-‘■i ~ ^il + ^jl + • • • ) 


(7) 

( 8 ) 


A’ = Af> + zw y;)”' + (-») 6mW 4'>"' 

+ [!!;>"' + 4'> I™] - f/'> [c™ + ?fi>] +,,,. (9) 


The symbols have the same meaning as in Eq. (29) of [13]. The ellipses denote further 
terms that do not contribute to the kernels i = 1, 2 of the current-current operators. The 
matrices and contain the UV counterterms from operator mixing. Compared to 
the calculation of the tree amplitudes, they have to be extended by the mixing with the 
penguin operators including the correspondent evanescent operators na. This implies, 
in particular, that the one-loop amplitudes must be computed including the 0{e) 
terms for all operators Qj, which mix with the current-current operators. Finally, one 
has to convolute the hard-scattering kernels with the LCDA, for which we adopt the 
conventional Gegenbauer expansion. 


3 The topological QCD penguin amplitude 


In this section we give the numerical results of the penguin amplitudes and and 
discuss the size and scale dependence of the new contribution. At LO, the penguin 
amplitude coefficients are given in the CMM basis by (W = 3, Cp = 4/3) 

®4LO “ [^3 + CpC^ + 16(05 + CpCo)]. (10) 


They are identical for p = u,c and independent of the LCDA. At NLO we have {L = 
In Sp = mp/ml, u = 1 — u) 


a 


4 ,NLO I Ci ,2 


as Cp 

4^ iV 


Oo- 


Oi 

2 W 


du 


2 2 

■^L + --G{sp-ie,u) 


</m 2 (m), ( 11 ) 


where we show only the terms from the current-current operators to illustrate the struc¬ 
ture of the result. Here 


G(s, u) 


2(12s-|-5 m — 3m In s) 2^(2s + u) ■C + 1 

-In- 

9m 3m 'C ~ 1 


( 12 ) 


is the one-loop penguin function with ^ = \/l — ^s/u. In practice, one then inserts the 
Gegenbauer expansion of 0Af2 (m) truncated at the second order to perform the integral. 
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Figure 2: The LO, NLO and NNLO values of a^inK) and a1{7rK) in the complex plane. 
The NNLO point includes a theoretical error estimate. 


The result is hnally expressed in terms of Wilson coefficients, quark masses and the 
Gegenbauer moments 

At NNLO the explicit expressions are involved, and we postpone the details to a future 
publication. Our hnal numerical predictions for the leading QCD penguin amplitudes 
are given as (for input parameters, see section 0]): 

<(7riF)/10-2 = -2.87 - [0.09 + 0.09i]vi + [0.49 - 1.32i]p, - [0.32 + O.Tlijp^ 


sp 


,0.434 J 


I [0.13]lo + [0.14 + 0.12*]hv - [0-01 - 0.05*]hp + [0.07]tw3 


+ f-1.94+°-32h 


<(7riF)/10-2 = -2.87 - [0.09 + 0.092]vi + [0.05 - 0.622]p, - [0.77 + O.SOzjp^ 


+ 


r r 


sp 


L 0 . 434 J 


{ |0.13]lo + |0.14 + 0.12i|Hv + 10.01 + 0.03i]„p + |0.07],„3} 


= (-3.34«;«) + (-1,0511>;S)». (14) 

In both equations the second line represents the spectator-scattering term, which for 
Tsp = 0.434 makes only a small contribution to a^. In the respective hrst lines, the number 
without brackets is the LO contribution, which has no imaginary part, the following 
two numbers are the vertex and penguin NLO terms, and the new two-loop NNLO 
contribution from the current-current operators Q 12 is fhe number labelled P 2 . We 
observe that the new correction is rather large. It amounts approximately to 40% (15%) 
of the imaginary (real) part of ol^nK), and 50% (25%) in the case of a\{TiK). Graphical 
representations of a\{'KK) are shown in Fig. [2] at LO, NLO and NNLO, where the NNLO 
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Figure 3: The dependence of the leading QCD penguin amplitudes a^{7iK) on the hard 
renormalization scale /i (form factor term only). Dashed, dashed-dotted and solid lines 
represent LO, NLO, and NNLO, respectively. 

point includes the theoretical error estimate^ The larger uncertainty of the imaginary 
part of al{7iK) is a consequence of the sensitivity to the charm-quark (pole) mass, for 
which we adopt the conservative range rric = 1.3 ± 0.2 GeV. 

The values ffT3|) . ffTTj) depend on the renormalization scale due to the truncation of 
the perturbative expansion and on hadronic parameters. The dependence on the renor¬ 
malization scale /i may be considered as a measure of the accuracy of the approximation 
at a given order in perturbation theory. This is shown in Fig. Elfor the form factor term 
contribution to aF^ijiK) up to NNLO. We observe a considerable stabilization of the scale 
dependence for the real part, but less for the imaginary part. This is explained by the fact 
that the imaginary part vanishes at LO. Hence only the hrst correction is now available 
and is, moreover, large. 

4 Phenomenology — direct CP asymmetries 

We now consider the new contribution to in the context of the full QCD penguin am¬ 
plitude and provide hrst results for some direct CP asymmetries. We defer the discussion 

^The LO and NLO numbers here as in the subsequent figure are not the same as (HSl), (HU) truncated 
to LO and NLO, because they employ Wilson coefficients C at LO and NLO, respectively. Moreover, 
consistent with previous LO and NLO calculations, they are computed in the operator basis as defined 
in [5S]. On the other hand, in (TT^ . (1141) NNLO Wilson coefficients in the CMM basis are used throughout. 
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of branching fractions to the more complete treatment including the two-loop matrix 
elements of the penguin operators Qs-q. 

We recall that in the QCD factorization approach the full QCD penguin amplitude 
consists of the parameters af, Og, and the penguin annihilation amplitude /Sf in the 
combination |5] 

al{MiM2) = ± + /3f(MiM2), (15) 

where the plus (minus) sign applies to the decays where Mi is a pseudoscalar (vector) 
meson. The hrst term, 04 (MiM 2 ), is the only leading-power contribution. Its real part 
is of order —0.03. The annihilation term is l/m{, suppressed and cannot be calculated in 
the factorization framework. Estimates based on the model dehned in [1] suggest that 
its modulus is also of order 0.03. While the magnitude of these two terms is largely inde¬ 
pendent of the spin of the hnal state mesons, the contribution from the power-suppressed 
scalar penguin amplitude r^^ag(MiM 2 ) is very small when M 2 is a vector meson, but 
larger than the leading-power amplitude for pseudoscalar M 2 . It interferes constructively 
for the PP hnal state, and destructively for VP. It follows from this brief discussion 
that the impact of a correction to 04 is always diluted in the full penguin amplitude. 
When M 2 = V, the computation of ascertains the short-distance contribution to the 
amplitude, and hence the direct CP asymmetry, but there is an uncertain annihilation 
contribution of similar size. When M 2 = P, there is another NNLO short-distance 
contribution from Og, which is difficult though not impossible to calculate, since it is 
power-suppressed. These features will be clearly seen in the analysis below. 

In the following we adopt the same values for the Standard Model, meson and form 
factor parameters as in Table 1 of PI with the exception of |14b/Vc6| = 0.085 ± 0.015, 
tb^ = 1.52 ps, 771^(2 GeV) = (90 ± 10) MeV, and = (190 ± 10) MeV. The decay 
constants, Gegenbauer moments and form factors involving kaons coincide with [9] , those 
involving K* mesons with 0 , except for {Q) = 0.39 ± 0.06. We note that the B- 
meson LGDA parameter is not important here, since the leading spectator-scattering 
contribution to the QGD penguin amplitude is colour-suppressed. 

In Fig. |4] we show the QGD penguin amplitude a'i{MiM 2 ) normalized to the sum 
of colour-allowed and colour-suppressed tree amplitude Q!i( 7 r 7 r) -|- Q! 2 ( 7 r 7 r)jl as was shown 
before in [5(19]. but now includes the NNLO computation for numerator and denominator. 
The NNLO result is represented by the dark point with error bars and corresponds to 
setting pyi = 0 in the annihilation model, which implies a small value of /^g. The nearly 
circular contours around this point show the variation of the theoretical prediction when 
the phase of the annihilation model is varied from 0 to 27r for hxed qa = 1,2,3 (inner 
to outer circles). The radius of the circle for = 1 leads to the estimate l/Sgl ~ 0.03 
mentioned above. The LO and NLO results are marked by diamonds without error bars. 
Despite the sizable NNLO correction to 04 as shown in Fig. [2], the difference between 
NNLO and NLO is small. This is a consequence of the “dilution” discussed above and a 
partial cancellation in the ratio of amplitudes. 

^For M 1 M 2 = nK, ttK*. For M 1 M 2 = pK, pK* we use the pp final state instead. Also, for pK* and 
p/3, only the longitudinal polarization amplitude is considered in the following. 




Figure 4: The QCD penguin amplitude &i{MiM 2 ) for the PP = ttK hnal state and its 
PV, VP, and VV relatives. The VV case refers to the longitudinal polarisation amplitude 
only. Shown are the theoretical predictions for the ratios a1{MiM2)/{ailnn) + a2(vr7r)) 
{pp instead of tttt in the lower row) and a comparison of extractions of the modnlus (rings) 
and phase (wedges) from data. Note there is no data for the CP asymmetry in the rate 
of the longitudinally polarized p~^K*~ hnal state. See text for further explanations. 
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The theoretical prediction can be compared to data, since the amplitnde ratio can be 
related to CP-averaged decay rates P and direct CP asymmetries. We discnss this for 
the PP case, from which the others can be inferred by obvions replacements. First, to 
very good approximation [5] 


al{7rK) 


Vub 

U 

1 - 

0 

1 

CKi(7r7r) -1- CK2(7r7r) 


Vcb 

Ik 

1 - 

In3 

1 

0 

1 _ 


( 16 ) 


which determines the grey rings aronnd the origin. The darker rings are dne to the exper¬ 
imental errors in the branching fractions and the lighter ones inclnde also the nncertainty 
of \VuhlVch\ (added linearly). To obtain the wedges we dehne ^|J to be the phase of the 
amplitnde ratio shown in the hgnre, and 


Q;i(7rJ’^) -|- a1{'KK) 
aii^rni) + CK2(vr7r) 


(17) 


We then hnd 


— simjj 


IvuTZ 
Re 7^ 


COsV' 


1 

2 sin 7 Re 7?. 




/tt Ptt+X- 

Ik \/2P^-^oP^-^o 


2lcp(7r+7C-). (18) 


In previons discnssions m the experimental error on the observables on the right-hand 
side and the error on 7 combined was large, so that it was jnstihed to assnme that TZ is real 
and to neglect the theoretical nncertainty on Re TZ, which mainly stems from the colonr- 
snppressed tree amplitnde a2{7m). This is no longer the case. The onter wedge now 
inclndes the theoretical nncertainty on TZ and 7 , which is added linearly to the pnrely 
experimental nncertainties (inner wedge). The middle wedge inclndes the nncertainty 
from 7 only. Note that flT^ has two solntions as shown in the hgnre, bnt the wedge that 
does not match the theoretical prediction is exclnded by T^^+k-< 1 . 

Since the NNLO correction to the amplitnde ratio tnrned ont to be small, we can 
reaffirm the conclnsions from [9] in the light of signihcantly improved data. The different 
magnitnde of the PP pengnin amplitnde vs. PV, VP and VV is clearly rehected in the 
data as predicted. There is reasonable qnantitative agreement as indicated by the error 
bars and the small onion-shaped regions corresponding to = 1. An annihilation con- 
tribntion of 0.02 to 0.03 seems to be reqnired, except for the longitndinal VV hnal states. 
The red sqnare in the hrst three plots of Fig. 0] corresponds to the theoretical prediction 
with Pa = 1 and the phase 0 a = —55° {PP), (pA = —45° {PV), 0a = —50° {VP) (see [1] 
for the dehnition of these qnantities), which is similar to the favonred parameter set S4 
of [5]. Only the CP asymmetry of the ttK hnal state now appears to reqnire a valne 
larger than Pa = 1 for a perfect ht. 

Moving to the observables themselves, we show in Table [T] the theoretical predic¬ 
tions for direct CP asymmetries, dehned as the rate asymmetry between B and B de¬ 
cays, together with the world average of experimental data (last colnmn), compiled from 
HFAG [2S]- We focns on the pengnin-dominated 6 —?■ s transitions of non-strange B 
mesons to vriP hnal states and their PV and VP relatives. We also show the CP asym¬ 
metry diherence 

5{tiK) = Acp{7i^K-) - Acp{7r+K-) (19) 
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1 1 cr+0.21+0.55 
-0.22-0.84 

n 88 +0.16 +1.31 
-0.17-0.91 

Q 4 q+0.09+1.09 
*^•^^-0.09-1.15 

-14111 

1 

* 

o 

A +0.25+0.60 

l-OU -0.26-0.47 

1 40+0.27+0.69 
-0.29-0.56 

Q qy+0.05+3.18 
-0.05-0.67 

-3.814.2 

1 

* 

o 

iq or+2.40+5.84 
iO.OO -2.70-5.86 

18 1fi+3-ll+ 7-79 

iO.iU-3 52-10,57 

1 ^ SI +3.01+69.35 
10.01 -2.83-15.39 

-6124 

tt+K*- 

1 1 1 o+2.00+ 9.75 
il.lO -2.15-10.62 

1 Q 70 +3-37 +10.54 

ly. 'n_3 80-11,42 

qq Qy+4.35+86.20 
ZO.Uf -4.05-20.64 

-2316 

71°K*° 

1 7 9q+3.33+ 7.59 
i 1 .ZO -3,00-12.57 

1 r 1 1 +2.93+12.34 
lO.ii -2 65-10.64 

9 1 (-+0-39+17-53 
-0.42-36.80 

-15113 


9 «Q+0.72+5.44 
-0.67-4.30 

1 r^+0.45+4.60 
-'-•^^-0.58-9.19 

7 9(-+1.21+12.78 
' •^'9-1.34-20.65 

17125 

A{TrK*) 

y -1 q +1.38 +3.38 
— ( --lo -1.28-5.35 

q ^ir +0.67 +9.48 
-0.59-4.95 

1 Qq +0.19 +4.32 
-‘■•'^'^-0.18-7.86 

-5145 

1 

o 

n qq+0.07+0.16 
-0.07-0.27 

n qq+0.04+0.19 
'^•^^-0.04-0.17 

Q qQ+0.06+2.28 

'^•'^'^-0.06-2.39 

-12 117 

p^K- 

1 Q 91 +3.42+13.95 
ly.Ol -3,61- 8.96 

A 17+0.75+19.26 
^•1 ' -0.80-19.52 

AO 7q+7.07+ 44.00 
^O. ro_7 62-137.77 

37111 

p+K- 

r 1 q+0.95+6.38 
-0.97-4.02 

1 rn+0-29+ 8.69 
i.OU-0 27-10.36 

9 C nq+4.43+25.40 
ZO.yo_4 gQ_75 63 

20 111 

O 

o 

o ^q+1.59+2.31 
-1.65-1.69 

q QQ+1.66+3.60 

O.yy-1 71-7 44 

14 49+0.08+19.49 
’9-+^ -0.08- 8.78 

6120 

S{pK) 

1/1 ly+2.80+7.98 
-2.96-5.39 

r ^7+0.96+10.86 
"I'll'-l.Ol -9.79 

1 7 80+3-15+19.51 
l'-°'9-3.01-62.44 

17116 

A(pR) 

q yr+1.62+4.78 
' ^-1.66-6.48 

in 84+1-98+11-67 
iU.O^_2 09- 9.09 

9 4 q+0.46+ 4.60 
^•+'9 -0.42-19.43 

-37137 


Table 1: Direct CP asymmetries in percent for the vriC, ttK*, and pK final states. The 
theoretical errors shown correspond to the uncertainties due to the CKM and hadronic 
parameters, respectively. The errors on the experimental values of 5 and A are computed 
from those of the individual observables appearing in fl 20 |) ignoring possible correlations. 


11 









and the asymmetry “sum rule' 


= Acp{7r+K-' 


^7r 4 


^ F' _n - 


( 20 ) 

The latter quantity is expected to be small [27] , since the leading CP-violating interference 
of QCD penguin and tree amplitudes cancels out in the sum. In order to focus on the 
effect of the new NNLO correction on the perturbatively calculable short-distance part 
of the CP asymmetry, the columns labelled “NLO” and “NNLO” give the respective 
results, when the long-distance, power-suppressed terms are set to zero. This means that 
we set /Sg to zero, as well as power-suppressed spectator-scattering terms. However, we 
keep the short-distance dominated, but power-suppressed scalar penguin contributions. 
The column labelled “NNLO-I-LD” adds the previously neglected terms back. The main 
effect is from weak annihilation, for which we adopt the S4-like scenario ( 54 ) marked by 
the red square in Fig. 01 

Focusing hrst on the “NLO” and “NNLO” results, we note that for the PP hnal 
states the change is minor, since, as discussed above, represents only part of the short- 
distance penguin amplitude. The situation is different for the ttK* hnal states where the 
Og contribution is small, and for the pK hnal states where due to the opposite sign of 
04 and Og a cancellation occurs. In these cases, we observe a large modihcation for the 
, 7r~^K*~ and the corresponding pK hnal states, for which the CP asymmetry arises 
predominantly from the imaginary part of d^/ai. These modihcations are a rehection of 
the sizable corrections seen in Fig. [2l The ehect is much less pronounced in the remaining 
modes, where the asymmetry is due to interference with (in case of p~K^) or 02 

(in case of p^K^), and the ehect of the NNLO correction cancels to a certain extent 

in the ratio of interfering amplitudes. Despite these large modihcations of some of the 
PV and VP modes’ asymmetries, the long-distance annihilation contribution is always 
more important numerically, and usually required to obtain a satisfactory description 
of the data. The modelling of the long-distance contribution also determines the hnal 
theoretical uncertainty, which can become very large. Given that the short-distance 
contribution is now known to NNLO and given the large amount of experimental data, it 
becomes imperative to better determine the annihilation amplitude, presumably through 
hts to data. 


5 Conclusion 

The computation of direct CP asymmetries in charmless B decays at next-to-next-to- 
leading order in QCD has been a long-standing issue. The long- and short-distance 
contributions can in principle be of the same order and a NNLO calculation is required 
to ascertain the perturbative part. In this paper we computed the two-loop contributions 
of the current-current operators 2 to the QCD penguin amplitude, which are expected 
to constitute the dominant contribution, at least to the imaginary part, which is required 
for observing CP violation. We hnd a sizable correction to the short-distance part of 
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the direct CP asymmetry, the effect of which is, however, tempered by power-suppressed 
short- and long-distance terms. Our preliminary conclusion is that the NNLO correction 
does not help resolving the nK CP asymmetry puzzle, nor does it render the poorly 
known annihilation terms redundant. The hnal analysis should, however, include the 
penguin operator matrix elements, as well as the one from the chromomagnetic operator 
considered in m. The corresponding calculations are in progress. 
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